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Abstract 

We confirm our previous limit on the free -space TiTi oscillation time T m i n — 
4.7 • W 31 y. The approach is verified by the example of exactly solvable model. 
We also discuss the all existing remarks. The approach with finite time interval 
includes potential model as an extremely simplified special case. It is shown that 
potential model is inapplicable to the problem in hand because for n-nucleus 
decay probability it gives W = 2. In particular, as we transparently show, this 
"decay law" (normalized by factor 2) is used by Dover, Gal and Richard for 
refutation of our results. 
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Recently, the calculations of nn-mixing in nuclear matter beyond potential model 
have been done[l]. It was shown that the nn- transitions in medium are suppressed 
only by factor 0.5, as a result of which 

r nn >T mm =T/V2 = A.7-10 31 y (1) 

where, T = 6.5 • 10 31 y[2]. This lower bound increases the previous one by 31 orders of 
magnitude. If this limit is true, it would rule out many existing GUTs. Certainly, this 
unusual result calls for comprehensive investigation. 

The aim of this paper is: l.To check the validity of the approach [1] by the example 
of the exactly solvable model. 2. To comment upon the main existing remarks. 3. To 
find out the conditions such that the restriction (1) is invalid. In doing so the reasons 
for disagreement with the previous results as well as the connection between S-matrix 
calculations and the approach with finite time interval are studied. 

1. Let us consider the nn transitions in nuclear matter (see Fig. 2 of Ref. [1]). The 
self-consistent neutron potential U n is included in the neutron wave function \n p >, 
V — (jPn/^ m + U n ,p n ). Block T is described by Hamiltonian //"(£) = (all n-medium 
interactions) — U n . For the whole process probability W(t) we get[l] 

W{t) = W a (t) + W b (t), W a {t) = e 2 t 2 (2) 

W\t) = -e 2 J* dt a fj dtpW n (t a , t p ), W n (r) = 2/m7f (r) (3) 

where W a ,W b are the contributions of Fig. 2a and 2b of Ref.[l] respectively, Wn(r) is 
the probability of n-nucleus decay in time r = t a — tp\ e = l/r nS . The matrix element 
TJ^ corresponds to n-nucleus decay: 

iT^tp) = f]H) fc < 0n p | /*" dh... f" 1 dt k H(t l )...H(t k ) | 0n p > (4) 

k=l Jt l3 Jt 

where | 0n p > is the state of medium containing the antineutron with 4-momenta 
p. In contradistinction to the matrix element of the whole process (see Fig.2 and 
Eqs.(8),(9) of Ref. [1]) this equation contains only Hamiltonian of n-medium interaction 
(annihilation) and represents a general expression for decay matrix. Two step process 
was reduced to decay of prepared state which follows the exponential decay law 

W n {r) = 2Imlf i {T) = W exp {r) = 1 - e~ rt (5) 
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where T ~ lOOMeV is the annihilation width of n-nucleus. It corresponds to all n- 
nucleus interactions followed by annihilation. However, the main contribution gives 
annihilation without rescattering of n[3], because a ann > 2<r sc . From (3) and (5) we 
get 



W b (t) = e 2 t 2 



1 1 1 

2 + Yt + TH i 



(e"»-l) 



(6) 



For our problem T = 6.5 • 10 31 y[2] and in Eq.(3) one can put W fi (r) = 1. Then 
W(t) = e 2 t 2 /2. If theoretical calculation of X^(t) is required, Eq.(5) will serve as a 
test. 

Let us consider the potential block T model, H = 5U = U, n — U n = const, where 
Un = ReUn — iT/2 is the n-nucleus optical potential. Direct calculation of T^(t) gives 



T? i (T)=i[l-exp(-iSUT)], 
W£°\t) = 2/mX*(r) = 2 fl - e" rr/2 cos(rRe5U) 



(7) 



Here W£ ot is the n-nucleus decay law obtained in the potential approach framework. 
Substituting Eq.(7) in Eq.(3), one obtains 



-e 2 t 2 



1 + 2Imi/(8Ut) 2 (exp(-i5Ut) + iSUt - 1) 



The whole process probability predicted by potential model is 

W pot = W a + W h vot = 2Imi(e/5U) 2 [l - iSUt - exp(-i5Ut)]. 
This expression coincides with the solution of Schrodinger equations: 



(8) 



(9) 



(id t + V 2 /2m - U n )n(x) = en(x), (id t + V 2 /2m - U n )n{x) = en(x). (10) 

Therefore: (i) The finite time approach was verified by the example of exactly solvable 
potential model. It is involved in Eqs.(2),(3) as a special case, (ii) The self-energy part 
£ = SU have been obtained dynamically, starting with singular diagrams (see Figs. 
3d,3e of Ref. [1]). Consequently, the potential model substantiation was given (from 
view point of singularities, other than physics! ). 

In the strong absorption region Tt ~| SUt \» 1 the potential approach does not 
correspond to the physics of decays. Really, let us return to the decay law (7). When 
r » 1/r ~ l(T 24 s, 



W? ot = 2. 
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The decay matrix T£(t) obtained by means of U n gives incorrect result for observable 
value W n . Taking into account that W is represented through W n we conclude that 
potential model is inapplicable to the problem under study. 
In the region | SUt \« 1 potential model has a sense: 

W?(T) = W exp {r) = Tr. (12) 

It remains to see the reason of enormous quantitative disagreement with our results. 
The formal answer is that the overall factor 2 in Eq.(7), compared to Eq.(5), leads 
to the full cancellation of the e 2 t 2 terms in Eq.(9). The strong result sensitivity was 
to be expected because when 5U = antineutron Green function G ~ 1/0 that is 
a consequence of zero momentum transfer in the e-vertex. q = is a peculiar point 
corresponding to transition between degenerate states. We will return to this question 
below. 

The standard difinition of decay law is 

W n (r) = 1- | U tt | 2 , (13) 

where U is evolution operator. Taking into account that unitarity condition UU + = 1 
must be fulfilled for U-matrix of any process we immediately get Eq.(3): W n {r) = 
2JmTj"(r). We would like to stress that it is sufficient to put W n {r) — 1. In any case 
W n < 1 and W > e 2 t 2 /2. 

As for potential model situation is obvious: for H = 5U = const from Eqs. (13) 
and (4) we immediately get Eq.(ll). 

Summing up we remind that potential model was used only for finite time approach 
verification. In the strong absorption region Ft » 1 the potential approach is invalid 
in principle and so, for realistic problem the model-independent method based on 
Eq.(5) has been proposed. Moreover, one can put W n = const < 1. In this case only 
unitarity condition is employed. 

2. It is very unpleasant to respond to hasty and erroneous criticism, however it 
should be made. 

Dover, Gal and Richard[4] assert that for realistic problem Eq.(7) must be used, 
other than (5). They calculate W(t) substituting Eq.(7) and certainly get W = W pot . 
On this basis they assert that potential model is correct. The fact that W^ ot = 2 
is ignored. By the other words Dover, Gal and Richard calculate W substituting for 
n-nucleus decay probability W n = 2. 

Dover at al. also claim that I am using the equation W = —2iTu. This expression 
is applied by them (see Eq.(2) of Ref.[4]) and not by us, because we know that W = 
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llmTn. At last I would like to remind that calculation presented in Ref. [4] have been 
done in our previous paper (see last but one paragraph of Ref.[l]), however conclusion 
was opposite. 

The strong absorption region Ft ~| SUt \» 1 is interesting for us. Sometimes it is 
proposed to eliminate Fig.2a (Eq.(8)) of Ref. [1]. In this case: (1) For exactly solvable 
model we have W pot = W* t = -e 2 t 2 < 0. For our calculations W = W b = -e 2 t 2 /2 < 0. 
(2) Free-space oscillations (p — > 0, H — > 0) are not reproduced. (3) At last, Fig.2a 
corresponds to the first term of the formal [/-matrix expansion and the iterative series 
for (10). 

Krivoruchenko's preprint [5] is completely erroneous. We will cite only two main 
points. (l)"...the quantity < TT + > has no probability meaning, along with the 
imaginary part of the T-matrix element —21m < T >". (See pg.6 of Ref. [5].) To our 
knowledge \Tji\ 2 has a probability meaning by definition of T-matrix. Krivoruchenko's 
assertion means that the whole theory of reactions and decays is incorrect. (2) The 
initial eq.(ll) of Ref. [5] must describe the nn-transition, annihilation. However, the 
l.h.s. of eq. (11) is free of n -nucleus interaction at all. The r.h.s. contains annihilation 
width T and coincides with potential model result. We would like also to get the result 
without calculation, but some difficulties emerge in reaching this goal. 

One additional comment is necessary regarding t-dependence of the whole process 
probability W(t). Eqs.(2),(3) have been obtained in the lowest order on e . The precise 
distribution W pr (t) which allows for the all orders on e is unknown. However, W is 
the first term of the expansion of W pr and we can restrict ourselves to a lowest order 
W pr = W, as it usually is for rare decays. 

The protons must be in very early stage of the decay process. Thus the realistic 
possibility is considered [6] that the proton has not yet entered the exponential stage 
of its decay but is, instead, subject to non-exponential behavior which is rigorously 
demanded by quantum theory for sufficiently early times. At first sight in accordance 
with (1) for nn-mixing in nuclear the similar picture should be expected. In fact the 
situation is more serious. We deal with two-step process. In attempting to calculate 
T in the framework of standard S-matrix theory we get T ~ 1/0. So there is no sense 
to speak about decay law exp(— Tt). It is necessary to calculate W{t) as it was done 
above. 

3. Let us try to render the result invalid. Assume for simplicity that in the nn- 
transition vertex (see Figs. 2a, 2b of Ref.[l]) the zero component of 4-momenta q is 
transferred. The antineutron Green function is 
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G(x - Xl ) = -ij ^^exp(ip n (x ~ xi) - i(t - ti)(p T 2 t /2m + U n + q ))6(t - h). (14) 



(2k) 

Instead of Eqs.(2),(3) we get 



W q (t) = W c (t) + W d (t) 



(15) 



W c (t) = 2Imi(e/q ) 2 [l - iq t - exp(-iq t)]. 



(16) 



W d (t) = 2e 2 / dt a I dtpe 
Jo Jo 



-irlmqo 



[ReT%(r) sm(TReq ) - 0.5W s (r) cos(rReq )} , (17) 



When g = 0, W q (t) = W(t). Let us consider Eq.(16). If go = SU then (16) coincides 
with W pot . In the region | SUt \» 1 

W c (t) -> W' c {t) = (2e 2 Iml/5U)t. (18) 

The S-matrix calculation [7] in the framework of the potential model (n appears in the 
state with U n ^ U n ) gives the same result 



W exp = T exp t = W' c (t). 



(19) 



Consequently, for nonsingular diagrams and big t finite time approach converts to the 
S-matrix theory: W c W' c = W exp . First term of W exp expansion is reproduced or, 
what is the same, T exp = dW exp /dt. This result is obvious and should be considered as 
a test. 

The previous lower bound r m j„ ~ ly is obtained from (18). Therefore, the energy 
gap g = SU 7^ leads to linear t-dependence (18) and, consequently, to the pro- 
cess suppression r min ~ ly. In our calculation g = and W(t) ~ t 2 . As a result 
W(T)/W exp (T) ~ TT ~ 10 62 . 

Eq.(17) analysis corroborates these conclusions. Really, let us take for simplicity 
ReT£ = Imq = 0. Eq. (17) becomes 



I r 2 1 

+e (^+r 2 ) 2 



W d (t) = e 2 [^(cos(g t) - 1) + 
e- rt (2g r sin(°g t) + (g 2 - T 2 ) co°s(g t)) + T 2 - g 2 



(20) 
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When Tt » 1, the linear t-dependence Wd ~ e 2 t is reproduced again. However, if 
go — ► then the principal contribution gives the first term and Wd ~ e 2 t 2 . 

The energy gap SU ^ is the reason for dramatic process suppression^! by the 
other words the restriction (1) is invalid when q ^ in the nn-transition vertex. 
SU 7^ can be realized by means of 3-tail (q ^ in the e- vertex), that is the other 
GUT; or effectively, in the framework of the potential approach. It is invalid in the 
region Tt » 1. 

The qualitative process picture is very simple. Both pre- and post-nn -transition 
spatial wave functions of the system coincide. Thereupon, n -nucleus decay takes 
place, which is described by T . The magnitude of T depends on system wave function. 
However, the result is T- independent, because when T ~ 10 31 2/, TT » 1 and W(T) = 
e 2 T 2 /2. When p — > 0, free-space oscillations are reproduced, which can be interpreted 
as a superposition of plane waves. In the slight absorption region Tt < (5 — 10) our 
results coincide with the potential model ones. 

In the next paper the following statements will be proved. l.The contribution of 
the corrections is negligible. 2. All the results are also true for any nuclear model. We 
will also draw special attention to the qualitative aspects of the problem. 

We do not see any effect which could change the restriction (1). 
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